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PHASE CONTROL OF SEISMIC DATA 

FIELD OF THE INVENTION 

5 [0001] This invention relates generally to the field of geophysical 
prospecting and, more particularly, to processing of seismic data. 
Specifically, the invention is a method for control of wavelet phase during 
seismic data processing. 

BACKGROUND OF THE INVENTION 

io [0002] Phase spectrum of the seismic wavelet is an important 
characteristic of seismic data. There are two major aspects of the seismic 
data phase control process. The first aspect is control of the seismic wavelet 
during intermediate processing steps. The second is phase control of the final 
seismic processing product, i.e., the seismic image. The first objective, 

is wavelet phase control during intermediate processing steps, is important 
because some data processing methods assume particular phase properties 
of input data. If the input data deviates from the particular assumption, the 
processing method will produce a related error. Thus, careful phase control at 
this stage insures that input phase specification errors are reduced and the 

20 efficiency of intermediate processing steps with respect to the wavelet phase 
treatment is enhanced. The second objective, phase control of the seismic 
image, is important because detailed interpretation of seismic data may 
depend on the wavelet phase, and interpreters usually prefer a seismic 
wavelet with a particular phase property such as zero phase. A zero-phase 

25 wavelet in the seismic image means that the data have the highest possible 
resolution for a given amplitude spectrum. Such data are usually easier to 
interpret. It should be noted that both of the above-described phase control 
aspects are related in the sense that better phase control during intermediate 
processing results in better quality data before final conversion to (for 

30 example) zero phase and hence to a better quality final (zero-phase) product. 
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[0003] For example, there are typically modules in the processing 
stream that assume that the data being processed have a minimum phase 
spectrum. Nevertheless, the data will typically have a phase spectrum that is 
not minimum phase. (A concise discussion of wavelet phase may be found in 
5 the Encyclopedic Dictionary of Exploration Geophysics by R.E. Sheriff, 4 th Ed. 
(2002) published by the Society of Exploration Geophysicists, in the definition 
of "phase characteristics." See also the definition of "wavelet" as well as other 
terminology used herein.) Thus, error might be introduced by processing data 
that is not minimum phase. 

10 [0004] One such module in a typical seismic data processing stream 
that assumes minimum phase property of the seismic wavelet is predictive 
deconvolution. Sheriff defines "deconvolution" as a "data processing 
technique applied to the seismic data for the purpose of improving the 
recognizability and resolution of reflected events." The purpose of performing 

15 predictive deconvolution is to remove a predictable part of the seismic data 
defined in terms of a prediction distance, thus attenuating periodic multiples 
and, as an option, compressing the wavelet. This technique is discussed in 
the literature, for example by Robinson and Treitel in Geophysical Signal 
Analysis, Prentice-Hall (1980). In discussing the method of predictive 

20 deconvolution on page 267, the authors state, "it depends on the deterministic 
hypothesis that the basic seismic waveform associated with each of these 
events is minimum delay." (See also page 29.) 

[0005] Another module commonly used in seismic data processing that 
assumes minimum phase, or delay, is the receiver consistent deconvolution 

25 method. The purpose of this method is to correct the wavelet shape for 
receiver-related effects. A related but more general technique is surface 
consistent deconvolution which provides decomposition of the seismic 
wavelet into source, receiver, offset and common midpoint terms. Application 
of only one term from the solution, the receiver term, is equivalent to the result 

30 from receiver-consistent deconvolution. 
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[0006] Gibson and Lamer pointed out the problem in 1984 in the 
context of vibroseis data processing (seismic data generated by a vibrator 
source). Predictive deconvolution is commonly applied to such data even 
though, as Gibson and Lamer state, "this process involves a minimum phase 

5 assumption." ("Predictive deconvolution and the zero-phase source," 
Geophysics 49, 379-397 (1984)) As the authors note, vibroseis data will be 
much closer to zero phase. They disclose a solution in the form of correcting 
the phase of the original vibroseis data toward minimum phase before the 
deconvolution, and they illustrate the effectiveness of this approach. They 

10 use a phase correction filter based on a correlated vibroseis signature, the 
Klauder wavelet. However, their method of designing a minimum phase filter 
relies on statistical inferences from the seismic data rather than a 
deterministic approach such as is disclosed in the present invention. 
Moreover, Gibson and Lamer do not account for phase distortions that can be 

15 caused by processing steps (e.g., frequency filters) that often precede 
predictive deconvolution in a typical processing stream. Hootrnan and Hart 
propose using the Gibson-Lamer aj)p/oach to compensate for phase 
differences arising when mixed sources (e.g., vibrpseis and dynamite) are 
used in the same survey. ('The realities of processing mixed-source seismic 

20 surveys," 68? Ann. Internal Mtg., Spe. Expl. Geophys., Expanded Abstracts, 
1436-1439(1998)). 

[0007] Connelly and Hart propose a method of making a phase 
correction post stack, after predictive deconvolution has been applied (pre- 
stack) to the data traces. (55 * Ann. Internal Mtg., Soc. Expl. Geophys., 

25 Expanded Abstracts, 491-495 (1985)) For the most part, this method 
addresses only the second aspect of the phase control process, i.e., phase of 
the final seismic image. Furthermore, their post-stack approach relies on the 
commutative property of linear shift-invariant operations. Therefore, this 
requires all processes prior to the stack to be linear shift-invariant. Also, like 

30 Gibson and Lamer, no attention is given to the problem of phase distortion 
from steps preceding predictive deconvolution. 
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[0008] Hart and Hootman point out that surface-consistent 
deconvolution is "subject to the same assumptions made when performing 
predictive deconvolution, which includes that there is a minimum-phase 
wavelet . . ." ("Achieving consistent and stable phase with mixed-source 
5 surveys," paper given at September 5, 2000 Technical Luncheon, available at 
web site http://cseq.ca/luncheons/200009/ ) They propose using the approach 
disclosed in the Connelly and Hart paper cited above to correct this problem. 

[0009] Martinez (U.S. Patent No. 4,646,274) and Galbraith (U.S. Patent 
No. 4,348,749) teach a method for using a vibrator source signature to correct 
10 for phase distortion introduced in standard processing of the data. They 
propose measuring the true ground force imparted to the earth by the vibrator 
to develop a phase correcting inverse filter. Like the method of Connelly and 
Hart, the filter is to be applied to the data after the processing. 

[0010] Although minimum phase is preferred for processing steps such 
is as predictive deconvolution, a zero-phase wavelet is preferred for the post- 
processing interpretation steps. (See Connelly and Hart, and also Sheriff and 
Geldart, Exploration Seismology, 2 nd Ed., Cambridge University Press, 181 
(1995). This is because zero-phase wavelets have the simplest shape and 
the highest peak for a given amplitude spectrum. Further, the peak occurs at 
20 the reflection time of the event. This alignment is important since the seismic 
wavelet generally broadens with increasing depth, with a zero-phase wavelet 
remaining symmetrical about the event time. 

[0011] Robinson and Treitel make an argument at page 251 that 
reverberation pulse-train waveforms will arrive at the receiver exhibiting 
25 minimum delay, "or at least approximately so." Nevertheless, the person 
skilled in the art of seismic processing will know that this assumption that the 
seismic data is minimum phase is often not accurate. The present invention 
provides a method of avoiding the error that can result from assuming the 
data have minimum phase spectrum. 
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[0012] Duren and Anderson disclose in U.S. Patent No. 5,384,752 the 
following: "Because of the shortcomings of the statistical data analysis 
methods, deterministic methods are generally recognized as superior 
methods of analyzing seismic data. Deterministic methods involve directly 

5 measuring the waveform of the source pulse . . . [and using] a source pulse 
function that is based on the actual pulse instead of an estimated wavelet." 
What is needed is a thorough, deterministic, pre-stack method of correcting 
wavelet phase spectrum to be consistent with assumptions in common data 
processing steps such as predictive deconvolution and also to shape the 

10 phase spectrum of the processed data for the following interpretation phase. 
The present invention satisfies these needs. 

SUMMARY OF THE INVENTION 

[0013] In one embodiment, the invention is a method for controlling 
wavelet phase spectrum in processing of any seismic data, which comprise 

15 using the measured or calculated source signature wavelet, after processing it 
through the same initial processing steps as are used on the data, to 
construct a shaping filter to convert the data to minimum phase, applying the 
filter to the data before processing steps that assume minimum phase, 
processing the minimum-phase seismic data, using the parallel-processed 

20 source signature to construct a second shaping filter to convert the data to 
zero phase, and applying the second filter to the processed seismic data. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0014] The patent or application file contains at least one drawing 

25 executed in color. Copies of this patent or patent application publication with 
color drawings will be provided by the Office upon request and payment of the 
necessary fee. 
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[0015] The present invention and its advantages will be better 
> understood by referring to the following detailed description and the attached 
drawings in which: / 

[0016] figure 1 is a flow chart showing the primary steps of one 
5 embodiment of the present inventive method; 

[0017] figure 2A shows actual seismic data after typical processing 
steps with no phase control; and 

[0018] figure. 2B shows the same data after processing by the present 
inventive method. 

10 [0019] The invention will be described in connection with its preferred 
embodiments. However, to the extent that the following detailed description is 
specific to a particular embodiment or a particular use of the invention, this is 
intended to be illustrative only, and is not to be construed as limiting the scope 
of the invention. On the contrary, it is intended to cover all alternatives, 

is modifications, and equivalents that may be included within the spirit and 
scope of the invention, as defihed by the appended claims. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0020] In the present invention, the source signature, S(t) % is used to 
derive a filter for converting the data to minimum phase before such steps as 

20 predictive deconvolution and receiver consistent deconvolution. The source 
signature is the output of a seismic source measured at or near the source 
before the seismic wave can be altered by the media it must traverse before 
arriving at a detector. In the case of a vibrator source, the signature may be 
measured by accelerometers mounted on the vibrator. In the case of an air 

25 gun, the signature may be measured by locating a hydrophone in the water 
directly beneath the air gun. In some cases, the source signature may be 
modeled with sufficient accuracy in lieu of making the more expensive direct 
measurement. 



WO 2005/019864 



PCT/US2004/015753 



-7- 

[0021] Figure 1 is a flow chart indicating the main steps in one 
embodiment of the present inventive method. At step 10 in figure 1, the 
source signature is processed through the same initial processing steps as 
the data from the seismic receivers. Those initial steps may include, for 
5 example, an anti-alias filter followed by bimating, or selecting every other 
sample of the digital data. The anti-alias filter is normally a high-cut filter 
applied before a sampling step such as the bimating step. The purpose of 
such a filter is to attenuate high frequencies that will "alias" (cause frequency 
ambiguity) upon sampling. The purpose of bimating is to reduce the data 

10 volume to be processed without seriously affecting the accuracy of results. 
Another example of an initial processing module that might come next is a 
band-pass filter to attenuate noise. Step 10 typically includes all processing 
modules that do not assume minimum phase and that come before the first 
module in the processing stream that does assume minimum phase. At a 

15 minimum, step 10 should not include any modules that do assume minimum 
phase if the full benefit of the present invention is sought. 

[0022] As stated, the (digitized) source signature wavelet is put through 
the same initial processing steps as the data. An exception would be a step 
that would not be applicable by its nature to the source signature. An 
20 example of that is the scaling and summation step that is taught by Duren in 
U.S. Patent Publication No. 2003/0048696 (Mar. 13, 2003) in order to remove 
the receiver ghosting interference from two-component ocean bottom cable 
data! 

[0023] At step 20 of figure 1, the source signature S(t) after step 10 is 
25 used to design and construct a filter that will convert the source signature to 
minimum phase. The mathematical basis for one such filter is as follows, 
where S mph (t) is the minimum phase signature equivalent to S(t). 

[0024] The Fourier transform of S(t), a complex quantity, may be written 
as 
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S(a>) = \S(<o)\exp{- j[Ps(a>)l , (1 ) 

where \S(o))\ is the signature's amplitude spectrum, and Ps(a>) denotes the 
signature's actual phase spectrum, <o is frequency and j = V^T . 

[0025] By definition, S mph (t) is a function with the same amplitude 
5 spectrum, \S(a>)\ t as S(t) but with a minimum phase spectrum, P mp h(o>). 

Pmph((o) and Ps(o)) are, of course, generally different. The problem of finding 
Pm P h((o) is known as "spectral factorization" which has been investigated in the 
art. (See J.F. Claerbout, Fundamentals of geophysical data processing, 
Blackwell Scientific Publications, 49 (1985)). Though the problem is not trivial, 
10 a number of spectral factorization methods have been proposed. For 
instance, the Kolmogoroff method establishes a relationship between the 
logarithm of the amplitude spectrum and the corresponding minimum-phase 
spectrum via Hilbert transform. Thus, P mp h(t) may be estimated from a 
logarithm of \S(g>)\ using Hilbert transform. 

15 [0026] The desired phase-shaping filter, F mph (t), can be written in 
frequency domain as follows 

(") = <*P{- j\P mP H (") " M (2) 

[0027] The amplitude characteristic of the filter is constant and equal to 
a unity. The person skilled in the art will know variations on and alternatives 
20 to the above-described method of designing a phase-control filter. All such 
methods are within the scope of the present invention. 

[0028] The filter constructed at step 20 is then applied in step 30 to the 
seismic data. This will shape the seismic wavelet to its minimum phase 
equivalent. In one embodiment of the present invention, the process is 
25 finished. The data is now suitable for processing by modules such as 
receiver-consistent deconvolution and predictive deconvolution that assume 
minimum phase. If for some reason the processing steps after step 30 
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assume some phase other than minimum phase, the filter designed in step 20 
would be constructed so that its application would bring the phase of the data 
to the appropriate phase. 

[0029] In another embodiment of the present invention, the previously 
5 described steps are followed by step 40 in which the remaining processing 
steps are performed on the seismic data and on the now minimum-phase 
source signature, i.e., S^ft), produced as a result of shaping by the same 
filter used on the data in step 30. At step 50, the source signature after step 
40 is used to design a second phase-control filter to shape the wavelet to 
10 whatever phase spectrum is desired for the data interpretation stage to follow, 
usually zero phase. A zero-phase wavelet is preferred because for a given 
amplitude spectrum, zero-phase wavelets have the simplest shape and the 
highest peak. Further, the peak occurs at the reflection time of the event. 
This alignment is important since the seismic wavelet generally broadens with 
is increasing depth, but a zero-phase wavelet will remain symmetrical about the 
event time. This second phase-shaping filter may be designed, in principle, in 
the same way as described above for the first phase-shaping filter. 

[0030] Let S'^it) be the source signature after step 40, i.e., after 

being subjected to the processing steps used on the seismic data such as 
predictive deconvolution, surface consistent decon volution, band pass 
filtering, etc. Thus, S' mph (t) and the output of the minimum-phasing filter of 

Eq.(2), S mph (t) % are generally not the same. One method of designing a 

phase-shaping filter to yield the zero-phase equivalent signature, 5'^ (*), may 

be described as follows. 

[0031] Analogous to Eq.(1), the complex spectrum of the S' mph (t) t 
S'mph (<*>). ma Y be easily calculated, e.g., using a Fourier transform, and written 
as 



20 



12 

25 

r 



S' api (<*>) = \S' mph {a>\exp{- (a>)J , 



(3) 
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where |S'^ A (^)| and F mph {a>) denote the amplitude and phase spectra, 
respectively. 

[0032] By definition, the zero-phase equivalent of S' mph (t) is a function 
S' 2pk (t) with the same amplitude spectrum, \s f mph , but zero phase 
5 spectrum, i.e., F 2ph = zero. The problem is trivial, and S' 2pfl (t) is completely 
defined by |S' m/ , A (a] . 

[0033] The desired shaping filter F^if) is a phase filter that can be 
written in spectral domain as follows: 

io F tp M-exp\j[F mp Ml (4) 

The amplitude characteristic of the filter is constant and equal to a unity. 

[0034] Just as with the minimum-phasing filter, persons skilled in the art 
is may know of variations or alternative ways to construct the zero-phasing filter. 
The filters from Eq.(2) and Eq.(4) have been derived in frequency domain. An 
alternative approach using least mean-square shaping filters may be used to 
produce finite length filters in the time domain (e.g., see Robinson and Treitel, 
170). 

20 [0035] At step 60, the second phase-control filter, e.g., the filter from 
Eq.(4), is applied to the data, which are now ready for (possibly) removal of 
phase dispersion due to attenuation and for subsequent processing and 
interpretation. A deterministic dispersion correction is generally required 
before the data are considered to have zero phase spectrum. This is typically 

25 accomplished by application of a filter called an inverse Q filter. To clarify the 
reference to subsequent processing, the processing flow addressed in this 
application or patent, including the appended claims, does not include such 
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possible steps as seismic imaging (migration) stack and post stack 
processing, for which the phase control preference is usually the same as for 
interpretation. 

[0036] The value of controlling the phase of the data to be consistent 
5 with the assumptions of the processing algorithms is common to all seismic 
data, P-wave or S-wave, however acquired. For example, seismic vibrator 
data is typically initially correlated with a reference signal, the pilot sweep 
signal that is used to control the vibrator. The result of such a correlation is to 
limit the phase spectrum to zero phase, or close thereto. This is not optimal 

10 for ensuing processing steps if those steps assume a different phase 
spectrum such as minimum phase. For a vibrator source, the source 
signature may be measured, for example with geophones or hydrophones 
located appropriately relative to the source, or with accelerometers mounted 
on the source, or the signature may be modeled. Similarly, data produced by 

15 impulsive sources such as explosives may not be minimum phase. 
Signatures for explosive sources may be measured or modeled. 

[0037] Figures 2A and 2B show actual examples of processed ocean- 
bottom-cable seismic data. The data in figure 2A were processed using a 
conventional processing stream, whereas figure 2B shows the same data 
20 after processing with the same processing stream but with the incorporation of 
the present inventive method. A clearer and higher resolution image is 
obtained using phase control as provided by the present invention. For 
example, the higher amplitude (red) bands in the range 1800 - 2000 ms are 
more sharply defined and less smeared. 

25 [0038] The foregoing description is directed to particular embodiments 
of the present invention for the purpose of illustrating it. It will be apparent to 
one skilled in the art that many modifications and variations to the 
embodiments described herein are possible. All such modifications and 
variations are intended to be within the scope of the present invention, as 

30 defined in the appended claims. 
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WE CLAIM: 

1. A method for use in seismic data processing to control the 
wavelet phase spectrum of seismic data collected using a seismic source with 
a known signature wavelet and at least one receiver, said method comprising 
the steps of: 

(a) using the source signature to design a first phase-control 
filter to limit wavelet phase spectrum consistent with assumptions inherent in 
at least one subsequent processing step, having first processed the source 
signature with any processing steps previously applied to the seismic data; 

(b) applying the first phase-control filter to the seismic data 
before said subsequent processing steps; and 

(c) performing all subsequent processing steps on the 
phase-controlled seismic data. 

2. The method of claim 1, wherein the first phase-control filter 
shapes the source signature to minimum phase. 

3. The method of claim 1, wherein said subsequent processing 
concludes with a second phase-control filtering of the seismic data, said 
second phase control filter being designed in steps comprising applying the 
first phase-control filter to the source signature from which it was designed, 
then processing the filtered source signature with said all subsequent 
processing steps up to the second phase-control filtering, and then using the 
resulting source signature to design a second phase-control filter to limit 
wavelet phase consistent with post-processing needs. 

4. The method of claim 3, wherein the second phase-control filter 
shapes the source signature to zero phase. 
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5. The method of claim 1 , wherein said at least one subsequent 
processing step comprises the steps of receiver-consistent deconvolution and 
predictive deconvolution. 
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Finish processing data except for 
final phase shaping; do same to 
source signature. 
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Apply second filter to seismic data. 



(57) Abstract: Method for controlling the phase spectrum of 
seismic data to match assumptions inherent in subsequent pro- 
cessing steps. The source signature, after processing with the 
same initial processing steps (10) used on the data, is used 
to design a phase control filter (20) that shapes the seismic 
data (30) to have a minimum phase spectrum or whatever other 
phase spectrum the subsequent processing algorithms may as- 
sume. The processed data is then filtered with a second phase- 
control filter (60), also designed (50) using the parallel-pro- 
cessed signature, to shape the data to zero phase or whatever 
other phase may be desired for interpretation of the data. 
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